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Observation of the Radiative Decay 07 
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Abstract 

We report the observation of the decay 07 with a statistical significance of 5Aa in 78.1 fb~"^ 

of data collected by the Belle experiment at the KEKB e+e~ collider. This is the first observation of 
a flavor-changing radiative decay of a charmed meson. The Cabibbo- and color-suppressed decays 
— >■ (/)7r°, (prj are also observed for the first time. We measure branching fractions B{D^ —>■ = 
[2mtoli (stat.) tofr (syst.)] x 10-^ B{D^ (pn^) = [8.01 ± 0.26 (stat.) ± 0.47 (syst.)] x IQ-'', and 
B{D° (prj) = [1.48 ± 0.47 (stat.) ± 0.09 (syst.)] x lO'^^. 

PACS numbers: 13.20.Fc,14.40.Lb 
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FIG. 1: Decay diagrams for the decay — s- (pj: (a) short-distance, (b) long-distance process. 

Flavor-changing radiative decays of the charmed meson system, D — > V'-f where is a 
vector meson, have not previously been observed. In the Standard Model, the short- distance 
contribution to these decays is negligible (the branching fraction from this contribution is 
predicted to be less than 10~^), and the long-distance contribution due to a vector meson 
coupling to a photon is expected to be dominant 0,01 • Examples of short- and long-distance 
processes for the decay D'^ — > 07 are shown in FigureJTl The branching fraction of this mode 
is expected to lie in the range (0.04 - 3.4) X 10-5 U 13; the current 90% confidence level 
(C.L.) upper limit from CLEO is B{D^ ^7) < 1.9 x lO"'* In the i?-meson system, 
where the short-distance contribution to i? ^ X^'y {X^ is p, u) decay can be used to measure 
the CKM matrix element Vtd, theoretical estimates of the long-distance contribution are also 
uncertain j^. Measurement of radiative D meson decays therefore provides an important 
constraint on the interpretation of -B — > X^'j results. 

In this paper, we present the first observation of the (pj decay. The analysis is 

based on 78.1 fb"^ of data collected at the T(45') resonance by the Belle detector 0] at 
the KEKB e~^e~ collider 0. KEKB is a pair of electron storage rings with asymmetric 
energies, 3.5 GeV for e"*" and 8 GeV for e~ , and a single interaction point. The Belle 
detector is a large-solid-angle general purpose spectrometer that consists of a three-layer 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- flight scintillation 
counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect mesons and to identify 
muons (KLM). 

According to a Monte Carlo (MC) study, the most important backgrounds to —>■ 07 
are the Cabibbo- and color-suppressed decays — ^ 07r° and (prj, which have not previously 
been observed. We therefore conduct a search for these decay modes as well. To reduce 
the combinatorial background, candidates are combined with to form D*~^ candi- 
dates. We calculate the difference in invariant mass, AM = Mqo ^+ — Mjjo , and require 

cand cand 

143.4 MeV/c2 < AM < 147.4 MeV/c^. 

Charged particle tracks are reconstructed in the SVD and CDC and required to be con- 
sistent with originating from the interaction region; the nearest approach of the trajectory 
to the collision point, which is determined run-by-run, is required to pass \dr\ < 0.5 cm 
and \dz\ < 1.5 cm, where dz is taken in the direction of the positron beam and dr is in the 
plane perpendicular to it. Particle identification (ID) likelihoods for the pion (£^) and kaon 
{Ck) hypotheses are determined from the ACC response, specific ionization {dE/dx) mea- 
surement in the CDC, and the time-of- flight measurement for each track. To identify kaons 
(pions), we apply a mode-dependent requirement on the likelihood ratio 71 = Ck I {Ck + C,,): 
a criterion oi IZ > 0.6 {IZ < 0.1) yields an efficiency of 86% (79%) for kaons (pions). 
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The rate at which pions (kaons) are misidentified as kaons (pions) under these criteria is 
8% (5%). We select candidates from K^K~ combinations in the invariant mass range 
1.01 GeV/c^ < Mx+K' < 1-03 GeV/c^, and use combinations from the mass sidebands 
i^sb), 0.99 GeV/c2 < Mk+k- < 1-00 GeV/c^ and 1.04 GeV/c^ < Mk+k- < 1-05 GeV/c^, 
to estimate the background under the (f) peak. 

Neutral vr mesons (tt^) are formed by pairing photons, each with energy > 50 MeV, 
and requiring a pair invariant mass within ±16MeV/c^ (~ 3a) of the nominal ir^ mass; 
the photon momenta are then recalculated with a mass constraint 0]. For rj ^ •yy and 

— >■ 07 reconstruction, photon candidates passing > 50 MeV are selected if no pairing 
with any other photon in the event {E^ > 20 MeV) yields an invariant mass that is consistent 
with the 7r° mass (the vr" veto). Otherwise, rj meson candidates are formed from pairs of 
selected photons and required to have an invariant mass within ±40MeV/c2 (~ 4a)of the 
T] mass. The photon momenta are then recalculated with an r] mass constrained fit 

We make the following requirements on laboratory momentum or energy: P^^o > 
750 MeV/c for (pn^, Pr, > 500 MeV/c for ^r], and E^ > 450 MeV for D° (pj. 

The D* momentum in the e~^e~ center-of-mass is required to satisfy P^* > 2.9 GeV/c for all 
modes; this criterion is optimised for the selection of 07 and is above the kinematic 

limit for D*^s produced in B decays at the T(4S) resonance. These criteria are determined 
using MC such that they do not introduce bias for the signal in data. 




M.,„(GeV/c2) 




FIG. 2: — > 07r^: (a) invariant mass distribution for data (points with errors), combinations 
of 0-sideband and 7r° (shaded histogram), and sideband + MC signal shape, scaled to the result 
of the fit described in the text (solid histogram); (b) fitted D yield in bins of cos^hei (points with 
errors) (^hei is the helicity angle) and the MC prediction (histogram). 
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FIG. 3: (pTj: (a) invariant mass distribution for data (points with errors) and combinations of 

(^-sideband and rj (histogram); (b) background-subtracted distribution, with a fit function described 
in the text (solid curve). 



The reconstructed M^^^o distribution shows a clear enhancement at the mass (Figure 
El^a)). The distribution from non-resonant — >■ K~^K~n^ decay also peaks at the D'^ 
mass. Its contribution is evaluated using 0-sideband candidates, constructed from 0s67r° 
combinations and found to comprise 3.1 ± 0.9% of the net yield. The signal yield is 
extracted by a fit to the sideband-subtracted M^^o distribution, assuming an exponential 
shape for the background and a signal shape developed by the Crystal Ball experiment j^. 
We obtain 1254 ± 39 events {x^/ndf = 25.6/38). To measure the helicity state, we form 
the helicity angle 6hei, the angle between the K'^ and 3-momenta in the rest frame of 
the (J) meson. Due to the conservation of angular momentum, the distribution in cos dhei is 
expected to be proportional to cos^ 9hei for (pn^ and tpr] decays but proportional to 

sin^ 9 hei (= 1 — cos'^ 9 hei ) for — > ^7. The distribution for 07r° candidates agrees well 
with the MC expectation (Figure |2tb)). When a 07r° or (prj decay is reconstructed as a 07 
candidate by missing one photon, the distribution of cos 9 hei is still close to cos^ 9 hei- The 
background to 07 from (j)T[^ / (prj decay is thus strongly suppressed by a requirement 

on 9hei- 

The M^^ and 0-sideband (M^^^^^) distributions are shown in Figure a). We extract 
the signal yield from the sideband-subtracted M^ri distribution shown in Figure Ol^b). The 

fit yields 31.1 ± 9.8 signal events, where Gaussian and first-order polynomial shapes are 
assumed for the signal and background, respectively. The significance of the signal, taken to 
be -^Z— 2 In (Lq / L^est) where Litest and Lq are the maximum likelihood values with the signal 
floated and fixed to zero, respectively, is AAa. 

For the radiative mode, we require | cos 9hei\ < 0.4. A clear peak is observed at the Z)° 
mass in the M^^ invariant mass distribution (Figure ID^a)). The signal yield is extracted 
using the binned maximum likelihood method. The shapes of signal and background from 
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— > 07r°, (pi] and — > 07r+7r° (feed-down) are obtained by MC simulation. The reliability 
of the simulation is checked by studying K^'j combinations with an invariant mass near 
the D° mass: the process -^s7 is forbidden, so any peaking in the signal region 

must be due to — > K^tt^ and K^r] events, similar to the feed-down background in this 
analysis. We observe no peaking. We find similar distributions in data and MC. The feed- 
down rates are normalized to the observed rates for the source modes, and the systematic 
uncertainty is due to the uncertainty on the source rates. The shape of the combinatorial 
background is estimated by fitting a first order polynominal to 0-sideband {4>sbl) candidates. 
The normalization of this contribution is allowed to float in the fit. The extracted yield is 
27.6^1'^^ (stat.) ^'^5^0 (syst.) events. The significance of the signal is 5.4o". 

With the signal region defined as [1.78 — 1.92] GeV/c^ and the requirement on cosOh^i 
released, the distribution in cos Ohei shows an excess of signal over background near cos Ohei ^ 
(Figure El^b)). It is clear that, without the requirement on helicity angle, the feed-down 
from — ^ 07r° and (f)r) is large. Figure Efc) shows the helicity angle distribution after 
background subtraction. It is consistent with the sin^ 9hei distribution expected for the 
D° — ^ 07 signal. 

To minimize systematic uncertainties, we measure the branching fraction as a ratio to 
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FIG. 4: 07: (a) invariant mass distribution for data (points with errors), the fit described 

in the text (solid histogram), the background component of the fit (dashed), ^vr*^ background (dark 
shading), and the sum of ^vr'', 0ry, and cpir'^TT^ backgrounds (light); (b) cos^hci distribution 

in the signal region, with the MC predictions: total (solid), total background (dark), and non- 
(j)7r^ background (light); (c) background-subtracted cos 6*1161 distribution and the MC prediction 
(histogram) . 
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— > K^K , where we have a signal of 21787 ± 226 events, and derive the branching 
fraction using the world average B{D^ -> K+R-) = (4.12 ± 0.14) x 10"^ [l^. Many of the 
systematic errors associated with tracking and particle ID are at least partially canceled in 
the ratio. 

The reconstruction efficiencies are estimated via MC simulation. Some differences in 
efficiency between data and MC have been noted, in particular for photon and 7r° recon- 
struction. For 77°, this is studied using the double ratio of two modes of the rj, 

gdata(27r°) ^ Naataiv ^ 3n^)/NMciv ^ 37r°) 
eMc(2vrO) Ndataiv 1i)/Nmc{t1 ^ 77) 

where -Re(7r° 77)1^^3^0 = Re{ri 77) {Re = Sdata/eMc) is assumed. The efficiency 
correction factors estimated from this study are found to be (97.9 ± 2.7)% for 7r° [E^ > 
50MeV and P^o > 750MeV/c), (96.2 ± 1.6)% for rj {E^ > 50MeV and > 500MeV/c) 
and (98.5 ± 1.9)% for a signal 7 with E^ > 450 MeV (assuming e^o = x e^). The overall 
detection efficiencies are summarized in Table lU The 7r° veto results in a low reconstruction 
efficiency for the (pt] mode. The branching fractions for the observed decays and their 
ratios to the reference mode are summarized in the Table IIIII 



TABLE I: Efficiencies for reconstructed modes [%] 



Mode {D^ ^) 


K+K- 






Reconstruction 


10.4 


6.48 


2.22 4.32 




±0.1 


±0.04 ±0.05 ±0.04 


B{(j) K+K-) 




49.20 


49.20 49.20 


B{Ti^/r] 27) 




98.80 


39.43 ~ 


Efficiency correction 




97.90 


96.20 98.50 


for / rj / ^ 








Total 


10.40 


3.08 


0.413 2.09 


TABLE IL Estimated fractional systematic errors [%] 






(1)7] — > (/>7 


Tracking etc. 


0.59 


0.59 


0.59 


Particle ID 


2.74 


2.74 


2.74 


AM 


0.98 


0.98 


0.98 


Mass of (j) 


0.94 


0.94 


0.94 


Efficiency correction 


2.75 


1.62 


1.94 


for -K^ /rj/^ 








Fitting & BG 


1.67 


2.01 


±2.46/-3.99 


B{D^ K+K-) 


3.40 


3.40 


3.40 


B{<t)^K+K~) 


1.42 


1.42 


1.42 


Bi-jT^/ri 27) 


0.03 


0.66 




MC statistics 


0.99 


2.76 


1.11 


Total 


5.88 


6.16 


±5.86/-6.65 
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TABLE III: Measured branching ratios B f /Bx+k- for decay modes (/) and branching fractions 
Bf = {B f /Bx+K-)^^K+K- ■ The first error is statistical and the second is systematic. 



Mode 


BflB^+K- 






Bf (xlO-^) 




ffi Q1 +1-70 +0.301 
[D.Oi 48 _o.36j X 


10 


-3 


[2.60 l°I?t°:l?]x 10-1 




[1.94 ±0.06 ±0.09] 


X 


10-1 


8.01 ±0.26 ±0.47 




[3.59 ± 1.14 ±0.18] 


X 


10-2 


1.48 ±0.47 ±0.09 



The systematic uncertainties are summarized in Table IHl The errors derived from the 
fitting process are estimated by varying the fit range, signal shape, and bin width in the 
fits; errors due to the backgrounds are estimated by varying background contributions within 
their uncertainties. The uncertainty in the acceptance due to the requirement on M^+k- for 

candidates is estimated by observing changes in the ratio Rfja^^^o = [A^^^ata/^A/cl^io^^^o 
while shifting the signal region by ±1 MeV/c^. Other systematic errors are estimated by ob- 
serving changes to the double ratio, Rdo^^t^o / Rdo^k+k- ■ The error due to the uncertainty 
in the tracking efficiency is estimated by changing the maximum impact parameter criteria 
from 0.5 cm to 0.3 cm for \dr\ and from 1.5 cm to 1.0 cm for \dz\. Similarly, the errors due 
to particle ID and the D* — mass difference are estimated by loosening the likelihood 
ratio cut to 7^ > 0.5 [TZ < 0.5) for kaon (pion) selection and by shifting the AM window by 

1 MeV/c^ on each side. Due to the difference in the momentum distribution between 
the — » 07r° and K^K~ modes, the uncertainty on the particle ID efficiency does 
not exactly cancel in the ratio. The uncertainties in the branching fractions of submodes 
are taken from the current world averages jl(i |. 

To summarize, we have observed for the first time a radiative decay of the D meson, 
in the mode — >■ 07. We also observe two other rare decays, — > 07r° and (prj, which 
are Cabibbo-suppressed and color-suppressed modes and constitute backgrounds for the 
radiative mode. The radiative D — » V'-y decays are expected to be dominated by long- 
distance contributions; however, the theoretical uncertainty on the rate is very large. The 
observed rate of — >■ 07 constitutes evidence for significant long-distance contributions, 
and it provides an anchor for the further development of non-perturbative QCD. 
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